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Abstract. We introduce and study the Hermitian matrix model with 
potential V s ,t(x) = x 2 /2 — stx / (1 — tx) , which enumerates the number of 
linear chord diagrams of fixed genus with specified numbers of backbones 
generated by s and chords generated by t. For the one-cut solution, 
the partition function, correlators and free energies are convergent for 
small t and all s as a perturbation of the Gaussian potential, which 
arises for st = 0. This perturbation is computed using the formalism 
of the topological recursion. The corresponding enumeration of chord 
diagrams gives at once the number of RNA complexes of a given topology 
as well as the number of cells in Riemann's moduli spaces for bordered 
surfaces. The free energies are computed here in principle for all genera 
and explicitly for genera less than four. 
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Introduction 

Consider a collection of b > 1 pairwise disjoint, oriented and labeled 
intervals lying in the real line RcC, each component of which is called a 
backbone. A chord diagram C on these backbones is comprised of a collection 
of n > semi-circles called chords lying in the upper half plane whose 
endpoints lie at distinct interior points of the backbones so that the resulting 
diagram is connected. 

This description of C with its chords in the upper half plane determines 
a corresponding fatgraph, namely, a graph in the usual sense of the term 
together with cyclic orderings on each collection of half-edges incident on a 
common vertex. The fatgraph C, in turn, determines an associated skinny 
surface £(C) with boundary by replacing each backbone with a rectangle to 
which are attached further semi-circular rectangles respecting orientation, 
one for each chord of C. In particular, S(C) contains C as a deformation 
retract; see Figure 1 for an example, where the chosen under /over crossings 
of the rectangles are immaterial to the fatgraph structure. 

chord diagram C as fatgraph skinny surface 2(C) 

Figure 1. From chord diagram to skinny surface. 

A chord diagram C therefore naturally determines an oriented surface 
S(C), which is connected since C is assumed to be connected. Thus, S(C) 
is characterized up to homeomorphism by its genus g > and number r > 1 
of boundary components, which we may equivalently associate to C itself, 
and the Euler characteristic of C or S(C) is given by b — n = 2 — 2g — r. 
This last relation can be used in particular to hnd the number of boundary 
components of the surface S(C) by inspection and then infer the genus. For 
example in Figure 1, we have b = 2 backbones, n = 4 chords and r = 4 
boundary components, which implies that the corresponding surface has 
genus g = 0. 

Let c 9t b(n) denote the number of isomorphism classes of chord diagrams 
of genus g with n chords on b labeled backbones whose generating functions 

(1) C gyb (z) = c gsb (n) z n , for g > 0, 

n>0 

are of our central interest here. We shall recursively calculate them using 
the topological recursion [15, 16, 21] of a Hermitian one-matrix model 

„ oo 

(2) Z= DH e - NtTV W = exp J2 N 2 ~ 2 9F g , 

J 9=0 
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where N denotes size of matrices, for a particular potential 



x 



2 



stx 



(3) 



V(x) 



2 



1-tx 



generalizing the Gaussian for st = 0. The free energy in genus g, which we 
compute, is given by 



where the constant terms reproduce the Gaussian free energies given by 
2g(2g-2) according to [24] in each genus g > 2 with appropriate modifications 
for g = or 1, and where E>2g denote Bernoulli numbers. The extra factor 
b\ arises because C g ^{n) counts chord diagrams with labeled backbones as 
opposed to unlabeled in the topological recursion, where a permutation of 
backbones must nevertheless preserve backbone orientations as we will see. 

As has long been known, cf. [6, 35, 27], Hermitian matrix models are 
well suited to various computations of combinatorial fatgraph generating 
functions which are determined by corresponding potentials. In section 2, 
we derive the specific matrix model potential V{x) = stx/{\ — tx) which 
encodes the solution to our particular combinatorial problem. This is a two- 
parameter perturbation of the Gaussian, and the parameters s and t in the 
expansion of the partition function Z s j are essentially generating parameters 
for the numbers of backbones and chords respectively. The entire solution of 
the model is encoded in the free energies F g in (2), which arise in the large 
N expansion of the matrix integral. Provided the 't Hooft parameter T is 
non-zero (and we shall ultimately take it to be T = 1 to restrict to the model 
of the form (2)), we show there is a unique continuous one-cut extension of 
the Gaussian which converges for small t and all s. It is this solution that we 
compute here and show reproduces more elementary enumerative results. 

Let us explain in more detail what the topological recursion is and how 
we can use it to determine the solution of the above matrix model. For a 
general Hermitian matrix model in its formal large ./V expansion, one can 
write down so-called loop equations, which are Ward identities or Schwinger- 
Dyson equations for certain multi- linear correlators wjf\pi, ... ,p n ) gener- 
alizing the resolvent. The leading order equation among those identities 
specifies a so-called spectral curve, i.e., an algebraic curve on which the re- 
solvent is well defined. It also turns out that all correlators Wn 9 \pi, . . . ,p n ) 
and loop equations they satisfy can be encoded entirely in terms of this spec- 
tral curve. These loop equations can be solved in a recursive way [15, 16, 21], 
and in this manner, free energies F g (for g > 2) are completely determined 

by correlators w[ 9 \p). 

We stress that this entire procedure requires just the knowledge of the 
spectral curve and not the details of the matrix model from which this 
curve was derived. An important achievement of Eynard and Orantin [21] 
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was to realize that one can use the recursive solution of loop equations to 
assign correlators Wn\pi, ■ ■ ■ ,p n ) and F g to an arbitrary algebraic curve, 
not necessarily of matrix model origin. On the other hand, it is guaranteed 
that F g computed for the spectral curve of a matrix model reproduce the free 
energies. For various interesting applications of this topological recursion to 
other problems related to moduli spaces of Riemann surfaces or algebraic 
curves, cf. [8, 9, 10, 11, 13, 24, 20, 21, 22, 31, 32, 33, 44]. 

It was already essentially known that in terms of the cut endpoints a < b, 
the eigenvalue density p(x) and the spectral curve y = M(x)^J{x - a){x - b), 
one can compute the leading free energies Fq and F± of (2), and we follow 
this procedure for our particular matrix model to derive the multiplicities 
c g,b( n ) of chord diagrams, for g = 0,1, in appendix A. The genus zero free 
energy can be found from a solution to the variational problem [12, 29], 
which leads to the following form of the effective potential 

T f 

(5) F = --U dxp{x)V{x)-\% 

where A is a Lagrange multiplier enforcing preservation of the number of 
eigenvalues to = J p(x)dx. The genus one free energy can be computed as 
in [3] (for generalization to more cuts see [14]): 

(6) F 1 = -1 log (M(a)M(b)(a - bf). 

Explicit results for genus zero and one in the model with potential (3) are 
presented in appendix A. These results reproduce known answers (up to 
three backbones, as given in (10), (12) and (18)), and more importantly 
provide generating functions for arbitrary numbers of chords and backbones. 
For example, we determine that generating functions of diagrams in genus 
zero and one, on four backbones, with an arbitrary number of chords take 
the form 

CM«) = ^<» + + 72,3 + 2448 ^ + „. 

(7) [ ' 
= 24^(715 + 7551, + 12456^ + 2096,3) = ^ + _ _ _ 



(1-4*) 



10 



To illustrate the combinatorial complexity of chord diagrams and appreciate 
precisely the generating functions which we derive, we show in appendix C 
by explicit enumeration that co4(3) = 72 and co4(4) = 2448 as predicted 
by (7). 

In this paper, apart from the side-note of determining the general Fq (by 
saddle-point methods) and F\ (due to Chekhov [14]) and analyzing them for 
our model in appendix A, our main result is a procedure based on topological 
recursion for calculating any desired one-cut free energy F g , for any g > 2, 
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on a sufficiently large computer. We implement it to find, for example, 

t\l-a) 2 

2 ~ 24(W 4 (1 - 5 - 4a + 3cx 2 ) 5 (l + 5 - 4a + 3a 2 ) 5 X 

x (^160(5 4 (1 - 3<r) 4 (l - of - 80<5 2 (1 - 3<r) 6 (l - a) 8 

+ 16(1 - 3ct) 8 (1 - a) 10 + 5 10 (-16 + 219ct - 462cr 2 + 252a 3 ) 

+ 105 6 (1 - 3o-) 2 (l - o-) 4 (-16 - 126ct - 423a 2 + 2286a 3 - 2862a 4 + 1134a 5 ) 

+ 5S 8 (1 - a) 2 (16 + 189cr - 2970a 2 + 9549a 3 - 11286a 4 + 4536a 5 )^) , 

where a = t(a + b)/2 and 5 = t(a — b)/2. This can be solved exactly in 
one of s,t and perturbatively in the other to extract specific C g< b(z) as we 
will explain. These methods of computation apply in all higher genera as 
well, and our similar but longer expression for F$ is given in (63). We 
again reproduce known results (up to three backbones, given in (12) and 
(18)) and then find generating functions for arbitrary numbers of chords 
and backbones. For instance as consequences of our computations, we find: 

144/ 

C 2) 4(z) = -, (38675 + 620648z + 2087808z 2 

(1 — 4z) 

+ 1569328/ + 134208/), 

48/ 

(l-4z) 16 

+ 1073841 1392z J + 4580024832z' 4 + 236239616z° 



C 3,4<» = „ "7.,M fi ( 53416125 + U94366915z + 6557325096/ 



144/ 

C 25 (z) = =3- (2543625 + 62424520z + 375044396/ 



2 



+ 671666053/ + 314761848/ + 18335696/), 
720/° 

C 3i5 (z) = ^(360380790 + 112750768652 + 95744892585/ 

(1 - 4z)~ 

+ 282797424880/ + 291167707410z 4 + 85497242928/ 
+ 3218434848/). 



Let us also comment on the profound import and wide scope of these 
results. The numbers c g j,(n) are of significance in computational biology 
because they describe the possible non-trivial topological types of complexes 
of several interacting RNA molecules as follows. Each backbone is identified 
with the sugar-phosphate backbone (hence the terminology) of a single RNA 
molecule oriented from its 5' to 3' end. If two nucleic acids comprising 
the b RNA molecules participate in a Watson-Crick basepair, then we add 
a corresponding chord taking care that chord endpoints in each backbone 



6 J. E. ANDERSEN, L. O. CHEKHOV, R. C. PENNER, C. M. REIDYS, P. SULKOWSKI 

occur in the correct order corresponding to the primary structure, i.e., the 
word in the four-letter alphabet of nucleic acids that determines the RNA 
molecule. In this way, a complex of interacting RNA molecules determines 
a chord diagram, and we demand that it is connected in order to guarantee 
an appropriate non-triviality of the interaction. Note in particular that 
nucleotides not participating in basepairs play no role in this model, i.e., 
there are no isolated vertices. 

The genus of the diagram determines the topological complexity of the 
interaction. This genus filtration has been profitably employed in a number 
of studies [34, 38, 42, 7, 23, 18] as well as in folding algorithms [41, 39] in 
the special case of a single RNA molecule (i.e., 6 = 1), which will be further 
discussed in section 1.2. Likewise, a folding algorithm for two interacting 
RNA molecules (i.e., 6 = 2) pertinent to antisense RNA for instance has 
been studied in [4]. 

At the same time, a simple transform of these numbers Cg^n) count a 
sub-class of chord diagrams called "shapes". These are discussed in the next 
section and give the number of cells in the ideal cell decomposition [37] of 
Riemann's moduli space for a surface of genus g > with 6 > 1 boundary 
components provided 2g + 6 > 2. 

The numbers computed here using the topological recursion are there- 
fore at once of significance in computational biology and in geometry thus 
representing a remarkable confluence of biology, mathematics and physics. 

This paper is organized as follows. There is background material in the 
next section including computations with the Gaussian potential itself for 
chord diagrams on at most three backbones. In section 2, we introduce 
our new model and determine the form of its potential (3). In order to 
solve this model, we first review some general relevant properties of matrix 
models in section 3 and introduce the formalism of topological recursion 
in section 4. Finally in section 5, we solve the matrix model (2) with the 
potential (3): we begin by determining its spectral curve and then compute 

the correlators Wn 9 \pi, ■ ■ ■ ,p n ) and free energies F g assigned to it through 
the topological recursion formalism. As the potential (3) depends on pa- 
rameters s (which generates numbers of backbones) and t (which generates 
numbers of chords), also the spectral curve, and then correlators and free 
energies depend on these two parameters as well. In appendix A, we present 
the computations of Fq and Fi, which are performed independently of the 
topological recursion, and in appendix B, we write down the (perturbative) 
solution of the equations determining cut endpoints in our model. In the 
final appendix C, we directly enumerate particular classes of chord diagrams 
just to confirm that indeed co,4(3) = 72 and co,4(4) = 2448. This enumera- 
tion also gives a sense of the simplest cases of objects we are counting here 
and can be read first if desired. 



TOPOLOGICAL RECURSION FOR CHORD DIAGRAMS 



7 



1. Background 

1.1. Chord diagrams, seeds and shapes. Two chords 7 and 7' in a chord 
diagram C with respective endpoints x,y and x',y' are said to be parallel 
if x, x', as well as y,y', lie in a common backbone with no chord endpoints 
in between, where x < x' and y' < y. Parallelism generates an equivalence 
relation whose equivalence classes are called stacks. 

Suppose the endpoints x,y of a chord 7 lie in a common backbone /3. If 
there are no chord endpoints between x, y, then 7 is called a pimple on /3, 
and if all chord endpoints in (3 lie between x and y, then 7 is called a rainbow 
on p. 

A seed is a chord diagram where every stack has cardinality one so that 
each pimple is a rainbow, and a seed is a shape provided every backbone has 
a rainbow. In particular for a chord diagram of genus zero, an innermost 
chord with both endpoints on a single backbone is necessarily a pimple, so 
for example, the only seeds of genus zero on one backbone are the empty 
diagram with no chords and the diagram with a single chord given by the 
rainbow, and only the latter is also a shape. 




seed but not shape shape 



Figure 2. Stacks, pimples, rainbows, shapes, seeds. 

Proposition 1.1. Other than the case g = 0, b = 1 (which was just dis- 
cussed), a seed of genus g on b backbones with n chords must satisfy the 
inequalities 2g + b — 1 < n < 6g — 6 + 5b, and these constraints are sharp. 

Proof. The lower bound on n follows from the fact that 2g + b—l < n+l — r 
according to Euler characteristic considerations together with the obvious 
constraint that r > 1. Conversely, a seed which saturates this lower bound 
necessarily has r = 1. 

If the skinny surface associated to a seed has more than one boundary 
component, then there must be a chord with different boundary components 
on its two sides; removing this chord decreases r by exactly one while pre- 
serving g again from Euler characteristic considerations. Define the "length" 
of a boundary component to be the number of chords it traverses counted 
with multiplicity. If there are ui boundary components of length £, then 2n = 



8 J. E. ANDERSEN, L. O. CHEKHOV, R. C. PENNER, C. M. REIDYS, P. SULKOWSKI 



J2e $- v i since each side of each chord is traversed exactly once in the bound- 
ary. It follows that In = 2{b + r + 2g-2) > v\ + 2v 2 + 3(r - v\ - u 2 ) > 3r-26 
since 2v\ + v 2 < 2b (except in the excluded case for which 2z/i + v 2 < 46), 
so 4(6 + g — 1) > r. There can thus be at most 4g + 46 — 5 such removals of 
chords to produce a seed with r = 1 proving the upper bound on n. □ 

Let Cg t b(n), s gi b(n) and t g ^{n), respectively, denote the number of iso- 
morphism classes of (connected) chord diagrams, seeds and shapes of genus 
g > with n > chords on 6 > 1 backbones. In each alphabetic case of 
X = C, S, T, we define corresponding generating functions 

X g,b( Z ) = E X 9,b( n ) 
(8) n>0 



and let 

(9) X g (z) = X g>1 (z). 

Note that whereas C g<b {z) is a formal power series, S 9tb (z) and T 9tb (z) are 
polynomials by Proposition 1.1. 

In particular, Cq{z) is the generating function for the Catalan numbers, 
i.e., co(n) = 00,1(77,) is the number of triangulations of a fixed polygon with 
n + 2 sides. They evidently satisfy co(n + 1) = J27=o c o(O c o( n — with basis 
co(n) = 1, which gives Cq{z) = 1 + z[Cq(z)] 2 , whence 



(io) cb W = 

Theorem 1.2. The generating functions for seeds and chord diagrams are 
related by 

(C Q {z)-r 



C g , b (z) = [C (z)]° S g , b ^ CQiz) 



S g,b{ z ) 



z + 1 
l + 2z 



b 

a 



g,b 



z(l + z) 
(l + 2z) 2 



Furthermore, the generating functions for seeds and shapes are related by 

(1+Z) b Tg, b ( Z )=Z b Sg, b (z). 

Proof. Writing simply Cq for Cq(z) and using Cq — 1 = zCq, we find 



= 1 % 1 1) = ^0 E(^o 2 ) 1 = £(^ 2 y 

2 -G l-(Co-l) 4 > 



As an argument of Sg tb , the jth term in the sum corresponds to inflating 
a single arc in a seed to a stack of cardinality j > 1 as well as inserting 
a genus zero diagram immediately preceding each of the resulting 2j chord 
endpoints. There is yet another factor Co arising from the insertion of a 
genus zero diagram following the last endpoint of the seed on each backbone 
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accounting for the further factor Cq. The resulting chord diagram is con- 
nected if and only if the seed itself is connected, and this proves the first 
formula. 

For the second formula, direct calculation shows that z = rfrPph inverts 



the expression u = 2-^Co(z) = ^2^j\~tz ' so ^ ne ^ rs ^ f° nmua reads 



S 9,b( u ) 



u(l + u) 
(l + 2n) 2 



-, -b 



a 



g,b 



u(l + u) 



Further direct computation substituting z = r^^^i into Cq(z) = 1 ^ 4z 

shows that Co( n~p^r ) = ^ipjj! an d the expression for S g fi follows. 

The third formula is truly elementary since a shape is by definition simply 
a seed together with a rainbow on each backbone. □ 

Seeds are useful insofar as chord diagrams as well as their generalizations 
allowing isolated vertices which arise in practice can be enumerated from 
them, cf. [41, 40] for instance. Shapes are important because of the following 
fundamental observation from [5, 37]: 

Theorem 1.3. Provided 2g + b > 2, the set of shapes of genus g > on 
b > 1 backbones with its natural partial order under inclusion is isomorphic 
to the ideal cell decomposition of Riemann's moduli space of a surface of 
genus g with b boundary components. 

Proof. Collapsing each backbone to a distinct vertex produces from a 
shape a fatgraph where each vertex has a loop connecting consecutive half- 
edges arising from the rainbow. The dual fatgraph of this fatgraph arises by 
interchanging the roles of vertices and boundary components. The collection 
of duals of shapes of genus g on b backbones gives precisely those fatgraphs 
that arise in the Penner-Strebel ideal cell decomposition of the associated 
Riemann moduli space, where the collapsed rainbows are dual to the tails 
discussed in [37]. Removal of chords corresponds to contraction of dual 
fatgraph edges. □ 

In light of Theorem 1.2, the explicit calculation of the generating function 
for chord diagrams we perform here thus gives as a consequence also the 
numbers of cells of given dimension in Riemann's moduli spaces of bordered 
surfaces. 

1.2. The special case of one backbone. For one backbone, there is the 
remarkable formula of Harer-Zagier [26] , which we discuss here, that arose in 
the calculation of the virtual Euler characteristics of Riemann moduli spaces 
for punctured surfaces and in principle solves for the various c g (n) = c 5i i(n). 
In effect, the Penner matrix model [35] directly computed the former without 
recourse to the latter. There is a large literature on the Harer-Zagier formula, 
cf. [28] and the references therein including several derivations of it. Here is 
this beautiful and striking formula: 
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Theorem 1.4. We have the identity 



Cg (n)N^-^ +1 _(l + z\ N 
>np< (2n-l)!! 

n>0 2g<n v ' 

Furthermore, the c g (n) satisfy the recursion 

{n + 1) c g (n) = 2(2n - 1) c g (n - 1) + (2n - l)(n - l)(2n - 3) c g -i(n - 2). 
This recursion directly translates into the ODE 

(11) z (i- Az )± Cg {z) + (l-2z)C g (z) = $ ff _i(z), 

where 

z" 2 <ViW = ^ 3 ^C 5 _!(z) + 24z 2 ^C 5 _!(z) + 27z£c g - 1 (z) + 3<Vi(s) 

with initial condition C s (0) = 0. From the ODE (11) and direct calculation 
follow several results: 

Corollary 1.5. For any g > 1, the generating function C g (z) is expressed 
as 



i • 

' 2 



(12) C g (z) 

(1 - Az) 3g - 

where P g (z) is an integral polynomial divisible by z 2g but no higher power 
and of degree at most (3g — 1) with P g (l/4) 7^ 0. 

The first few P g {z) are given as follows: 

Pi(z) = z 2 , 

P 2 (z) = 21z 4 (z + l) 

P 3 (z) = llz 6 (l58z 2 + 558z + 135) , 

P A {z) = U3z 8 (2339 z 3 + 18378 z 2 + 13689 z + 1575) , 

P 5 (z) = 88179z 10 (l354 z A + 18908 z 3 + 28764 z 2 + 9660 z + 675) . 

It is natural to speculate that the polynomials P g (z) themselves solve an enu- 
merative problem and that there may be a purely combinatorial topological 
proof of Theorem 1.4 based on a construction that in some way "inflates" 
these structures using genus zero diagrams, to wit 

C g (z) = P g {z){Vl^~z) 1 -^ 
Co(z) 



-C (z 

Pgiz^Coiz)) 69 - 1 (l + (C (z) - 1) + (C (z) - l) 2 + • • • ) 
P g (z)(C (z)) 6g ^ (l + z(C (z)) 2 + z 2 {C {z)f + • • • \*" ' 



63-1 
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Another calculational consequence of the ODE (11) is an explicit recursion 
for the numbers of seeds: 

Corollary 1.6. We have S±(z) = z 2 + 3 z 3 + 3 z 4 + z 5 , and for any g > 1, 
the following five-term recursion for the number of seeds: 

6 

(13) (n+l)s g (n) = 9h{n) s g -i(n - h), where 

h=2 

g 2 (n) = (2n - l)(n - l)(2n - 3) 
53 (n) = 2(4n-ll)(n-l)(2n-3) 
g 4 (ra) = (24n 3 - 180n 2 + 411n - 285) 
g 5 (n) = (16n 3 - 156n 2 + 464n - 414) 
g e (n) = 4(n - 2)(n - 4)(n - 6). 

The first few seed polynomials are thus given by: 

S 2 (z) = 21 z 4 + 210 z 5 + 840 z 6 + 1785 z 7 + 2205 z 8 + 1596 z 9 
+630 z 10 + 105 z 11 

S- 3 (z) = 1485 z 6 + 28413 z 7 + 225379 z 8 + 1006929 z 9 + 2862783 z 10 
+5496645 z 11 + 7325703 z 12 + 6812289 z 13 + 4348344 z 14 
+1819818 z 15 + 450450 z 16 + 50050 z 17 

S 4 (z) = 225225 z 8 + 6687252 z 9 + 83809440 z 10 + 605300410 z 11 
+2867032740 z 12 + 9542753220 z 13 + 23243924704 z 14 
+42438380985 z 15 + 58817592405 z 16 + 62093957640 z 17 
+49660516620 z 18 + 29612963952 z 19 + 12768025270 z 20 
+3763479720 z 21 + 678978300 z 22 + 56581525 z 23 

S 5 (z) = 59520825 z 10 + 2458871415 z 11 + 43395443091 z 12 + 447114000333 z 13 
+3067654998408 z 14 + 15065930976096 z 15 + 55383949540920 z 16 
+156912145081692 z 17 + 349331909457531 z 18 + 618798832452801 z 19 
+878354703383157 z 20 + 1001499812704755 z 21 + 915317216039226 z 22 
+665807616672198 z 23 + 380384810323518 z 24 + 166997907962886 z 25 
+54384464384250 z 26 + 12376076963250 z 27 + 1756856351250 z 28 
+117123756750 z 29 



Remark 1.7. [25] With further work starting from this recursion, one can 
show that the number of seeds is log concave, i.e., for 2g + 1 < n < 6g — 2, 
we have 

s g {n) 2 > s g (n + 1) s g (n - 1). 

In particular, it follows that the sequence s g (2g), s g (2g + 1), . . . , s g (6g — 1) 
is unimodal. 
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1.3. The Gaussian resolvents for at most 3 backbones. Provided the 
number b of backbones is less than four, the loop equations [2, 19] for the 
Gaussian potential permit the computation of the c g ^(n) in terms of c g (n) 
as we discuss here. This in part explains an important aspect of the com- 
putations of the previous section. See also [30]. 

For any s > 1, the s-resolvent or s-point function of the potential V{x) is 
defined to be 

(14) u s ( Xl , ...,*.) = N°~ 2 <tr — L_ . . . t r ^^) conn , 

where ($} denotes the part of the correlator J DH e~ NtrV ^ H ^(H) 
arising via Wick's Theorem from connected diagrams computed as a formal 
power series using tr -^jj = X)fc>o ^ r ~1ett> where the integral is over the 
N x N Hermitian matrices with their Haar measure 

(15) DH = f Ai<j<j<jv dReHij A dlmHij ^ A f Af =1 dff^ 

and with the partition function Z = J DH e~ NtrV ( H ) . 

For a polynomial potential V(x) = Ylk=l ^k xk f° r convenience (our meth- 
ods apply more generally whenever the potential has rational derivative), 
we find that 

g>0 m,...,n s >0 x l ' ' ' Xs 

e n (-tk) £k ^i. jns ({hA,... 



= E^ 23 ^ S 9) (X 1 ,...,X S ), 

where Uni,...,n s ({£i, ^3, • • • >^d}) denotes the number of isomorphism classes 
of connected chord diagrams of genus j on s labeled oriented backbones, 
with respective numbers n±, . . . , n s of incident half-chords, plus unlabeled 
oriented backbones with k incident half-chords, for k > 1. 

Now specializing to the Gaussian potential V(x) = \x 2 and letting ui s (x) = 
lo s (x, . . . , x) for convenience, we find 

u,(x) = E N~ 29 oj^(x) = x~ s E N~ 2 ° C g , s (x- 2 ) 

in the earlier notation since C g ^{z~ 2 ) = z s uj l f\z). The first Schwinger- 
Dyson equation or Ward identity with equal arguments in this case reads 
uj\ (x) — xloi(x) + 1 + N~ 2 uJ2{x) = 0, from which we solve for the 2-point 
function with equal arguments in terms of the 1-point function: 

(16) u 2 (x) = iV 2 {x^i(x) -u)x{x) - 1}. 
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Likewise, xu]%(x) = 2u}\(x)u!2(x) + ^u"(x) + N 2 ijj^{x) gives the 3-point 
function with equal arguments 



(17) 



1 



w 3 (x) = N {xlo 2 (x) - -oj'l{x) - 2ui(x)uj2{x)} 



in terms of u\{x) and its derivatives. 

Remark 1.8. One cannot continue this recursion to express higher resol- 
vents uj s (x), for s > 4, with equal arguments in terms of 0J\(x) alone because 
higher Schwinger-Dyson equations include products of traces (for example, 
ua(x) involves (tr -, — 3-— tr —^n) ) which cannot be evaluated in this 

*v / \ (x—H)- 3 x—H I conn' 

manner. We wonder if this is perhaps related to the behavior [36] of posets 
of chord diagrams: in the case of genus zero for s backbones with s = 1, 2, 3, 
the geometric realization of this poset is homeomorphic to a sphere; for 
s = 4, it is homeomorphic to a simply connected manifold which is not a 
sphere; and for s > 4, it is an increasingly singular non-manifold. 

Theorem 1.9. For any g > 0, we have 

PgM 



(18) 
where 



Cg,2(x) 



(1 - Axf9+ 2 



and Cg^(x) 



P 9 ;s(x) 
(1 - Axfs+ 5 ~¥ 



Pg,2( x ) = x lp g+i( x ) ~ Ph{x)P g+ i-h{x) and 



h=l 



P, 



gM x ) 



3+1 

^Pg+i^x) ~ 2 Ph(x)Pg+i-h,2(a 

h=l 



P 



9+1 



1 + (42 + 60g)x + (96 + 336c/ + 288g 2 )x z 



P' g+1 (x)(l-4x) 
P"(x)(l-4xf 



5x + (20 + 48g)x' 



2x' 



Proof. It follows from (16) that u)2(z) is given by 



N 



2 9J 9) I 



g>o 



Extracting the coefficient of N 2g , we conclude that 

9+1 



ZUJ 



(9+1), 



i.e., 



h=0 



9+1 



Z 0J 



z- 2 C g ,2(z~ 2 ) = zz^Cg+iiz- 2 ) -Y,z- 2 C h {z- 2 )C g+l _ h {z~ 2 ). 

h=0 
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Setting x = z~ 2 , we find 



9+1 

C g>2 (x) = x~ l C g+ i{x) - C h (x)C g+ i- h (x) 

h=0 

a 

= x- l C g+x {x) - ]T C h {x)C g+X - h {x) - 2C (x)C g+1 (x), 
1 1=1 



and rewriting this using (12) gives 

1 1-VT^X~\ Pg+l(x) 



X X 



{1-Axf9 +?J - L 2 



1 9 



(1 - 4x) 3 9+ 2 ^ 

(1 - 4x)39+ 2 



YPh{x)P g+ i- h {x) 



as was claimed. 

Starting from (17), the analogous but slightly more involved computation 
produces the asserted expression for C g ${x). □ 

For later comparison to our one-cut solution, we here record the first 
several such polynomials on two or three backbones in low genus: 



Pq,2{x) = x, 

P li2 (x) = x 3 (20x + 21), 

P 2>2 (x) = x 5 (l696a; 2 + 6096a; + 1485) , 

P 3i2 (x) = x 7 (330560a; 3 + 2614896a; 2 + 1954116a; + 225225) , 

P , 3 (aO = 2x 2 (3 + 4x), 

Pi l3 (x) = 12x 4 (45 + 207a; + 68a; 2 ), 

P 2 , 3 (x) = 6x 6 (15015 + 137934a; + 197646a; 2 + 27592a; 3 ), 

P 3>3 (x) = 8a; 8 (3132675 + 46335375x + 143262162x 2 
+98362965a; 3 + 7468348a; 4 ), 
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2. The potential V(s,t,x) 

If C is a disjoint union of chord diagrams, then let 

6(C) denote its number of backbones, 
n(C) denote its number of chords, 
r(C) denote its number of boundary cycles, 
Aut{C) denote its automorphism group 

permuting oriented backbones, and 
g(C) denote the sum of genera of its components. 

Furthermore, set = 2 n / 2 -k n2 / 2 and for any tuple t>i,t>2, • • • ,vk of non- 
negative integers, define 

P Vl ,..., VK (s,t,N) = — — L- f e- tr H ' 2 ' 2 l[(s tr (tff)*)"* DH 

for the integral as before over the N x N Hermitian matrices with respect 
to the Haar measure DH given in (15). 

Lemma 2.1. For any tuple vi, . . . ,vk, parameters s,t and natural number 
N, we have 

_ AT(C) s b(C) t 2n(C) 

Pv u ...,v K (s,t,N) = ^2 #Aut{c) , 

where the sum is over all isomorphism classes of disjoint unions of chord 
diagrams with backbones having k incident half- chords, for k = 1, . . . , K. 

Proof. The special case s = t = 1 is entirely analogous to that of Theorem 
2.1 of [35], where the replacement of ^ tr H k there by tr H k here corresponds 
to replacing fatgraphs there by chord diagrams here, i.e., distinguishing one 
sector at each vertex of the fatgraph to represent the location of the backbone 
in the chord diagram kills the cyclic permutations about the vertices. The 
general case then follows easily since b(C) = J2k v k and 2n(C) = J2k k y k- d 

A standard computation as on p. 49 of [35] thus gives 
Theorem 2.2. For Z(s, t,N) = -±f e~ N tT V ( S ^ H ) DH with potential 

(19) V(s, t, H) = hi 2 - = \H 2 -s £(tff)* 

we have 

N 2-2g(C) s b(C) t 2n(C) 



log Z(s,t,N)=J2' 



#Aut(C) 

where the sum is over all (connected) chord diagrams C. 



16 J. E. ANDERSEN, L. O. CHEKHOV, R. C. PENNER, C. M. REIDYS, P. SULKOWSKI 



An amusing point, which however does not help with computations, is 
that the derivative of the logarithmic potential from [35] 

■£-(x + log(l -x))=l + = —^—r 

dx\ J x — I x — 1 

agrees with our potential here up to the Gaussian factor for s = t = 1. 
Let us also emphasize that the partition function in this section enumerates 
chord diagrams with unlabeled backbones in contrast to the considerations 
in section 1.2. (Vertices from insertions are labeled whereas vertices from 
the potential are not in Wick's Theorem.) 



3. Matrix models and spectral curves: generalities 

As we sketched in the introduction, the main idea of this paper is to 
solve the problem of counting chord diagrams (or RNA complexes or cells in 
Riemann's moduli space) using matrix model techniques and the topological 
recursion. We present this general formalism in this and the next section. 
We will apply it in section 5 to find the free energy expansion of the one-cut 
Hermitian matrix model with the potential given in (19). 

Let us consider a general Hermitian matrix model 

/oo 
DHe-^ y W = exp ]r^ 2 F„ 

9=0 

with the measure DH introduced in (15). This model is slightly more general 
than (2) as we replaced the explicit size of matrices N in the exponents by 
1 /h. In this section, we present how to find a spectral curve of such a model 
in the large N limit, with the so-called 't Hooft coupling 

(21) T = hN = const 

kept fixed. In the next section, we will explain how, from the knowledge of 
the spectral curve, one can use the topological recursion to determine free 
energies F g . For a matrix model such as (20), the spectral curve and in 
consequence the F g depend on T. To solve the combinatorial enumeration 
problem in (2), we can simply set T = 1. 

The spectral curve arises from the leading order loop equation of a matrix 
model, and it can be characterized as a curve on which the resolvent is well- 
defined. The (all-order) resolvent is defined as a correlator 

, oo 

(22) u^x) = ^Tr— ^) = E h 29 4 9 Hx)- 

9=0 

The leading order term in this expansion u;[ (x) is often referred to simply 
as the resolvent. It is simply related to the density of eigenvalues p(x) which 
becomes continuous in large N limit and has the following asymptotics for 
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large x 

(23) J°\x) = I / '^U' ~ l - 

to J X — X' X 

where 

(24) t = / /o(x)da;. 



In this large limit, the eigenvalues are distributed along compact intervals 
called cuts, along which in particular the above integrals are performed. 
Moreover, an inverse relation determines the density of eigenvalues as a 
discontinuity 

(25) p{x) = —(ujf\x-ie)-uf\x + ie 

2m V 

in the resolvent along these cuts. The consistency of the above relations 
requires that 

(26) uf\x - ie) + J°\x + ie) = ^M. 

The conditions (23) and (26) determine the resolvent oo^f\x) completely, 
and it is often possible to guess its form by requiring that these conditions 
are met. Equivalently, one can find the resolvent using the so-called Migdal 
formula. In case there is a single cut with endpoints denoted a and b, this 
formula reads 



(27) J°\x)- 1 I dz V '( z )y/( x - a )( x - b ) 



2T Jc 2m x - z yj(z - a)(z - b) ' 

where the contour C encircles the cut. Moreover, expanding the resulting ex- 
pression around infinity and imposing the condition (23) provides a relation 
between 't Hooft coupling T and the cut endpoints. This one-cut resolvent 
will be sufficient for our computations in this paper, however, it is not hard 
to generalize the above formula to the case with many cuts having endpoints 
aj and bi\ in this case, one only needs to replace the radicals in numerator 
and denominator of the integrand respectively by rii^ — — b{) and 

Ui(z - <h)(z - bi). 

It is also convenient to encode the information about the resolvent u/°)(x) 
in a new variable 

(28) y = lv'[x) - Tuj(x) = ^(ujf\x - ie) - J°\x + ie)), 
which turns out to be related to x by a polynomial equation 

(29) A(x,y) = 0. 

This equation therefore defines an algebraic curve, which is called a spectral 
curve of a matrix model. As we will review below in the context of the 
topological recursion, one in fact needs to provide a parametric form of 
this equation, i.e., to represent it in terms of two functions x = x(p) and 
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V = y{p) °f a parameter p living on the Riemann surface such that the 
equation A(x(p),y(p)) = is satisfied identically. For the one-cut solution, 
the functions x{p) and y{p) are rational, while equation (28) takes the form 

(30) y = M(x)yJ(x-a)(x-b), 

where M(x) is a rational function. 



Example 3.1. Let us consider the Gaussian model with the potential 

l, 

2' 



V(x) = \x 2 . The Migdal formula (27) results in 



(31) ^ 0) (x) = ± f (x-sJ{x-a){x-bj) = 

, , a + b (a — b) 2 „ . 9 . 

and from the first line, it is clear that (26) is satisfied. The expansion in the 
second line together with the condition (23) then imply that 

(33) a = -2y/T, b = 2VT. 

Also, introducing y as in (28), we find the following equation for the spectral 
curve 

(34) y 2 = l( x -a)(x-b). 

Comparing with (30), we see that M(x) = ^ in this case, and a parametric 
form of this algebraic curve can be given for example as follows 



(35) 



4. Topological recursion 

The topological recursion introduced in [15, 16, 21] can be used to as- 
sign to an algebraic curve A(x, y) = a series of multi-linear differentials 
wjf\pi, . . . ,p n ) and functions F g called free energies. Free energies F g are 
symplectic invariants, i.e., they are invariant with respect to transformations 
which preserve the two-form dx A dy up to a sign. The multi-differentials 
transform in an appropriate way under symplectic transformations, and they 
depend on a particular parametrization x = x(p),y = y(p) of the algebraic 

curve. The differentials Wn\pi, ■ ■ ■ ,p n ) are defined recursively by equa- 
tions which, by definition, coincide with loop equations of a Hermitian ma- 
trix model [2]; as shown in [19], all details of these equations depend only 
on the form of its underlying algebraic curve A(x, y) = 0. Thus, applying 
the topological recursion to the curve which is a spectral curve of some par- 
ticular Hermitian matrix model such as (20), reproduces free energies F g 



TOPOLOGICAL RECURSION FOR CHORD DIAGRAMS 



19 



and correlators wjf\p\, . . . ,p n ). From the matrix model viewpoint, these 
correlators are defined as 



V -fa^)> 



x(pi)-M) \x(p k ) - M) / conn dx{p\) ■ ■ ■ dx(j) n ) ' 

which generalizes the resolvent introduced in (22), and in particular, (x) : 

ui^f\x)dx. For h = iV _1 , the functional part of Wn also agrees with ujn^ 
introduced in (1.3). 

The recursion relations for multilinear differentials Wn^ are defined in 
terms of the following quantities. Firstly, these relations are expressed in 
terms of residues around the branch points 

(36) dx(p*) = 0. 

For a point p in the neighborhood of a branch point p* there is a unique 
conjugate point p, defined as 

(37) x{p) = x(p) . 

Two other important ingredients are the differential 1-form 

(38) uj(p) = (y(p) - y(p))dx(p) 

and the Bergman kernel B(p,q). The Bergman kernel B(p,q) is defined as 
the unique meromorphic differential two-form with exactly one pole, which 
is a double pole at p = q with no residue, and with vanishing integral over 
Af-cycles § A B(p, q) = (in a canonical basis of cycles (Aj, B 1 )). For curves 
of genus zero as arise from the one-cut solution, the Bergman kernel takes 
the form 

(39) B M = J** 

VP - Q) 

A closely related quantity is a 1-form 



dE q {p) = l -J q B{^p). 



defined in a neighborhood of a branch point q* in terms of which we may 
express the recursion kernel 

(40) *<->^. 

With the above ingredients, we can now present the topological recur- 
sion. This recursion determines higher-degree meromorphic differentials 
Wn 9 \pi, . . . ,p n ) from those of lower degree. The initial data for the re- 
cursion are the following one- and two-point correlators of genus zero 

(41) Wi°\p) = 0, W^( Pl , P2 ) = B( Pl , P2 ), 
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with Wn = 0, for g < 0, by definition. For two sets of indices J and 
N = {1, . . . , n}, let us now respectively denote pj = {pi}i£j, and pjy = 
{pi, . . . ,p n }. The multi-linear correlators are then defined recursively as 



wii^PN) =J2 Re S K(q,p)(w^ 2 1) (q,q,p N ) + 



(42) 



m=0 JcN 



where Xjctv denotes a sum over all subsets J of N. In particular, the lowest 
order correlators determined by the recursion are 

Wl{p) = £ fos JT(g,p)W£(<z,g) 

and then [1, 17] 

W$(p, Pl ,p 2 ) = £ ResK((Z,p)(VF 2 °((Z,pi)^ 2 (g,P2) + W?(9,pi)W£(g,p2 

Finally, we can define the free energies F g in genus g. For g > 2, they 
come from the corresponding W^: 

(43) F g = ^_^Res0(g)^( g ), 

where 0(g) = J 9 y(p)dx(p). The expressions Fo and F\, given in (5) and (6), 
can be found independently as discussed, e.g., in [14, 21]. 

Example 4.1. Let us apply the formalism of the topological recursion to 
Example 3.1, i.e., the algebraic curve y 2 = j(x—a)(x—b) with a parametriza- 
tion given in (35). We find two branch points p* = ±1 and a global ex- 
pression for the conjugate point (valid around both of these branch points) 
p = p~ l . As the Gaussian curve represents the one-cut solution on the genus 
zero spectral curve, the relevant Bergman kernel is simply given by (39), and 
the recursion kernel (40) takes the form 

(44) K(q,p) = {a _ b)2{q2 _ l){p _ q){qp _ iy 
Computing the correlators, we find for instance 

w<%)~- 16p3dp 



(a-b) 2 (p 2 - 1) 4: 
^ 8d Pl dp 2 dp 3 2 2 

Wg I (Pl.P2.P3) =— ( 7T2- + (P2 + 1) (P3 + 1) 

(45) \ a — °) 

-(Pl-1)- 2 (P2-1)- 2 (P3-1)' 2 ) 

r ( 2 ) 86016dpp 7 + 3p 9 +p n 



ma 
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Finally, the genus-two free energy following from the definition (43) reads 

< 46 > F * - -TH^W 

Substituting a,b = ^f2y/T according to (33) we find 
(«) «— 

Computing higher order free energies, we find [24] well-known expressions 
for the Euler characteristics [26, 35] of moduli spaces of closed Riemann 
surfaces of genus g generalizing the above result for genus two, namely, 

(48) F g - — ^ l — 



2g(2g-2) T^-r 



5. Solution of the model 



In this section, we derive a formal large N expansion and find free energies 
F g of the matrix model (20) 



(49) Z = I DH e-^ v ^ = exp ]T h 2 ^ 2 F g , 

9=0 



with the potential given in (19), i.e., V(H) = ^H 2 — j^j- The potential 
in this model is a function of parameter s which generates the number of 
backbones and t which generates the number of chords. Thus, the F g are 
also functions of these two parameters. Expanding F g as a series in s and 
t we obtain multiplicities labeled by the genus, number of backbones and 
number of chords, which provide the solution of our model. In fact, free 
energies depend also on 't Hooft coupling T = HN, however, this dependence 
is rather simple, and all the information we seek can be obtained by setting 
T = 1. To solve this matrix model, we follow the steps explained in previous 
sections: we find a parametrization of the spectral curve and then apply 
the formalism of the topological recursion to derive associated multi-linear 
correlators wi 9 \p±, . . . ,p n ) and free energies F g . Finally, we present explicit 
expressions for free energies for small values of g and show that they indeed 
encode the c g fi(n) discussed in the introduction, i.e., the numbers of chord 
diagrams for various small genera and numbers of backbones. 

5.1. Spectral curve. We begin by determining a specific spectral curve. 
Let us note that the potential (19) is in fact not bounded below, and there 
is a singularity as x — > Nonetheless, this singularity does not pose a 
problem in the formal expansion of the free energy as we will see. Also for 
this reason, we treat the parameters s and t in the potential (19) as generat- 
ing parameters for powers of x, and so we effectively consider the behavior 
of eigenvalues only in the vicinity of the perturbed Gaussian minimum even 
though there are three critical points of the potential, i.e., we are studying 
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the one-cut solution of the model. We can therefore determine the one-cut 

,(°)/ 



semiclassical resolvent u}\ (x) using the Migdal integral (27) leading to 
(50) 



Coi , s x(xt — I) 2 — st J (x — a)(x — b) 

W j (x) = — ; " tts h — " X 

1 v ; 2T(xt - l) 2 2T 



st 2 (4 + (xt 2 - 3t) (a + b) - 2xt + 2abt 2 ) 



2{xt- l) 2 ((at- l)(bt- 1)) 



3/2 



(52) 



where a and b denote cut endpoints. Expanding this expression for large x 
we find that it has a form lo^\x) = c\ + C2X~ 1 + 0(x~ 2 ). However, imposing 
the asymptotics given in (23) requires that c\ = and C2 = 1. These two 
equations take the following form 

f = a + b+ ( ^+^) 
I ((«t-l)(W-l)) 
[ ' I , r r r _, ,vi , 4s((2-^^)(at+bt-2)+2abt 2 -3t(a+b)+i) 

iDi —{a-t>) h -372 . 

[ ((at-l)(W-l)) 

In principle, one could now solve (51) for a and b and substitute these values 
back into (50) to get the exact formula for the resolvent, but we know of no 
closed-form solution. However, one can determine perturbative expansion 
of a and b in s and t, and such a representation is sufficient at least as far as 
perturbative expansion of free energies is concerned. For example, we find 
that the perturbative expansions in s (with exact dependence on t) start as 

b = b(s,t) = 2TV 2 + st ^\-; 2 /2 +... 

and an analogous expansion in t (with exact dependence on s) can be de- 
termined. 

It is convenient to introduce half the sum and difference of the cut end- 
points a, b 

a + b a — b 

(53 S=^^, D = . 

v ' 2 ' 2 

Taking the square of the first equation in (51) as well as determining the 

denominator ((at — l)(bt — 1)) from both equations in (51) leads to the 
following pair of equations 

. A . j As 2 t\2,D 2 - 4T) 6 = D 2 {D 2 - AT) 2 (AD 2 - t 2 {3D 2 - AT) 2 ) 3 , 

( 54 ) \ c _ D 2 -4T 

I ° — t{3D 2 -4T) ■ 

The first equation is a polynomial equation for D which can be solved pertur- 
batively to any accuracy in s or t, and then the second equation immediately 
determines S. In appendix B, we present a perturbative expansion of S and 
D in s with exact dependence on t; from this one can easily find higher order 
terms in the expansion (52). 
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Theorem 5.1. Provided the 't Hooft parameter T / 0, there is a unique 
continuous extension of the Gaussian to a solution of (54)- In fact, the ex- 
tension is smooth in t and converges for all s and small t, and it furthermore 
satisfies a(s,—t) = —b(s,t). 

Proof. Consider the polynomial q S) t(S,D) given by the difference of the 
left- from right-hand side of the first equation (54). Direct computation 
shows that 9qs ^' D) ^ = 512 £> 15 (8T - 3D 4 )(AT - D 4 ). Provided this is 
non- vanishing, the Implicit Function Theorem applies to provide the unique 
smooth and convergent one-cut solution. Meanwhile for the Gaussian with 
t = 0, we have cut endpoints ±2y / T, so S = and D = 2y/T. Thus, 
T 7^ implies D ^ 0, and D 4 = 16T further implies that the factors 
(8t — 3D 4 ) = —8T and (4T — D 4 ) = — 12T are also non- zero as required. 

The asserted symmetry a(s,t) = —b(s, —t) is equivalent to the statement 
that S is an odd and D an even function of t. Meanwhile, the first equation in 

(54) involves only even powers of t so the smooth solution from the Implicit 
Function Theorem must be an even function of t. Insofar as T is independent 
of t, it is even as well, so the second equation in (54) displays S as an odd 
function of t. □ 

In what follows, we will work with general expressions for the resolvent 
and ultimately derive exact expressions for F g in terms of S and D. Sub- 
stituting the above expansions, we will then find perturbative expansions of 
F g in s and t. Let us therefore proceed to find the spectral curve in terms 
of parameters a and b or equivalently S and D. This curve is to be ex- 
pressed in terms of variables x and y defined as in (28). In particular, only 
expressions involving square roots \JJx — a)(x — 6), given in the second line 
of (50), will enter the formula for y. Moreover, note that the second line of 

(50) involves one non-rational expression ((at — l)(bt— 1)) , which however 
can be replaced by a rational function using the first equation in (51). In 
consequence, the equation for the spectral curve takes form 

1 / (tx - 1 + ( Q+fe ) t ) 2 + 7 

( 55 ) V = g \/(x-a)(x-b) ^ _ A , 

so that comparing with (30) we find 

The entire dependence in (55) on x and y is given explicitly, and 7 is a 
rational function of a, b and t given below. If we square this equation we 
obtain an algebraic equation for the spectral curve. In fact, notice that all 
factors of t in this equation multiply a and 6, and if we rescale also x and y 
by t, then the entire dependence on t can be encoded in rescaled variables 

(57) £ = xt, rj = yt, a = at, (3 = bt, a = St, 5 = Dt. 
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In particular in terms of variables introduced above, the expression for 7 
takes the form 

(a + /3)(a 2 + P 2 + 14(a + /3 - a/3) - 16) 
7 = 



(58) 



16(a + /3 - 2) 
(7(4 - 45 2 -7a + 3a 2 ) 



4(a-l) 

Finally in these rescaled variables, the remaining equation reads 

(59) (£ - 1)V = i(f - «)(£ - /3) [(f - 1 + ^) 2 + 7 

This is our desired expression for the spectral curve. In order to apply 
the topological recursion, we must further find a parametric form of this 
equation. To this end, note that the square root of the above equation, 
or equivalently the equation (55), is a modification of the Gaussian curve 
y = \\f{x — a)(x — b) given in (34) by a rational function in x. The essence 
of the parametrization of the Gaussian curve given in (35) was to take care 
of this square root and devise a function x = x(p) such that j(x — a)(x — b) is 
a complete square of a rational function in p, which is then identified with y. 
We can therefore borrow this aspect of the Gaussian parametrization (35) 
and just modify the form of y by multiplying it by this additional rational 
function of x in the present case. This leads to the following parametric 
form of the spectral curve 



(60) 




te-l+^) a +7 _ 5(1 n \ fe-l+f) 2 +7 

(FT? ~ Hp p ) (5-i) 2 ' 



(M- 



With this parametrization of the spectral curve, we are prepared to compute 
the free energies using the topological recursion. 



5.2. Free energies. We can now determine free energies F g , for genus 
g > 2, of the spectral curve given in (59) as parametrized by (60). This 
curve has genus zero, so the Bergman kernel takes the form (39). As the 
parametrization of the £ variable is the same as in the Gaussian case (35), we 
find the same branch points p = ±1 and globally conjugate points p = p~ l 
as given in Example 3.1. Nonetheless, all other quantities which depend also 
on r] have of course a much more complicated form than in the Gaussian 
case. In particular, the recursion kernel takes the form 

((g- 1 + q)5 + 2a - 2) 2 
v-v "y*,*; _ i)(p_ g )(gp_ 1)52 4 7 + (( g -i + ^ + 3(J _2) 2 ' 

With these ingredients, we can compute multi-linear correlators Wn^ ■ They 
are given by rather complicated and not very enlightening formulas, which 



/K1 \ lf( r, 



2q 3 



TOPOLOGICAL RECURSION FOR CHORD DIAGRAMS 25 

are too lengthy to write down here. However, just to present the simplest 
example, let us give the formula for : 



W$°\pi,p2,P3,) =2dpidp 2 dp, 



a - 1 
a z 



(62) 



a -5-1 



(Pi + 1) 2 (P2 + 1) 2 (?>3 + 1) 2 (1 + 5-4a + 3a 2 ) 
a + S-1 

' (pi ~ 1) 2 CP2 - 1) 2 (P3 - 1)2(1 - J - 4a + 3a 2 )- 



At this stage viz. Theorem 5.1, we can also confirm that all our results 
are consistent with the Gaussian model discussed in Examples 3.1 and 4.1. 
The potential (19) reduces to the Gaussian potential when st = 0. In 
both cases s = and t = 0, the system (51) simplifies so that the first 
equation implies a = —b, and then the second equation leads to the result 
(33) also in agreement with the expansion (52). For a = —b, we also find that 
7 = 0, and the curve (59) reduces to the Gaussian curve (34). Furthermore, 
the recursion kernel (61) and all correlators and (62) in particular reduce 
respectively to (44) and (45). 

Finally with Mathematica [43], we compute genus two and genus three 
free energies. As these are the free energies of the curve (60), which involves 
variables x and y rescaled by t, let us denote them by F 2 and F%. We find 
they are given by the following exact expressions 



240<5 4 (1 - 5 - 4cr + 3o- 2 ) 5 (l + 5 - 4a + 3o" 2 ) 5 
x (l60(5 4 (l - 3<r) 4 (l - af - 805 2 (1 - 3cr) 6 (l - a) 8 

+ 16(1 - 3a) 8 (l - a) 10 + <5 10 (-16 + 219a - 462a 2 + 252a 3 ) 

+ 10<5 6 (1 - 3a) 2 (l - a) 4 (-16 - 126a - 423a 2 + 2286a 3 - 2862a 4 + 1134a 5 ) 

+ 55 8 (l - a) 2 (16 + 189a - 2970a 2 + 9549a 3 - 11286a 4 + 4536a 5 )) 
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6 ~ 20165 8 (1 - 5 - 4a + 3a 2 ) 10 (l + 5 - 4a + 3a 2 ) 10 
x (l07520o 8 (l - 3a) 8 (l - a) 12 - 61440<5 6 (1 - 3a) 10 (l - a) 14 

+ 2304(W 4 (1 - 3a) 12 (l - a) 16 - 512(W 2 (1 - 3a) 14 (l - a) 18 + 512(1 - 3a) 16 (-l + a) 20 

+ <5 20 (512 - 48603cr + 471708a 2 - 1645251a 3 + 2623320a 4 - 1957284a 5 + 555660a 6 ) 

+ <5 18 (1 - a) 2 (-5120 - 639867a + 20457090a 2 - 177935466a 3 + 697802346a 4 

- 1439492067a 5 + 1623708828a 6 - 948940272a 7 + 225042300a 8 ) 

+ 252<5 10 (1 - 3a) 6 (l - a) 10 (-512 + 2970a + 70785a 2 - 3780a 3 - 1030050a 4 

+ 2252880a 5 - 1077300a 6 - 1682100a 7 + 2417850a 8 - 1129950a 9 + 178605a 10 ) 

+ 84<5 12 (1 - 3a) 4 (1 - a) 8 (1280 - 17523a + 256203a 2 + 5113206a 3 - 22143348a 4 + 2780703a 5 

+ 124222086a 6 - 273630879a 7 + 265968360a 8 - 126660105a 9 + 24111675a 10 ) 

+ 95 w (l - a) 4 (2560 + 156387a + 605388a 2 - 68079487a 3 + 696887576a 4 - 3203644311a 5 

+ 8142105636a 6 - 12207659661a 7 + 10771521264a 8 - 5182076088a 9 + 1050197400a 10 ) 

+ 3£ 14 (1 - 3a) 2 (l - a) 6 (-20480 - 6699a - 23380056a 2 + 166860351a 3 

+ 520822512a 4 - 6709872897a 5 + 22111131168a 6 - 36421137675a 7 

+ 33054947856a 8 - 15849836352a 9 + 3150592200a 10 )) 

The free energies F g in our model are closely related to these results. Indeed, 
the original curve is written in terms of x and y while the above free energies 
are computed for the curve in rescaled variables (57). In general, rescaling 
a variable x or y by t results in rescaling of the genus g free energy by 
t 2 9~ 2 . Because we rescale both x and y, this means that the free energies are 
rescaled by t 2(29_2) . The genus two and three free energies, with dependence 
on t reintroduced, are thus given as 

(63) F 2 = t A F 2 , F 3 = t 8 F 3 . 

Note that the entire dependence on t arises upon substituting a and 5 ac- 
cording to (57). 

The expressions (63) are our desired results, and all we need do is sub- 
stitute the expressions for a and 5 or in fact a and b given in principle by 
solving the equations (51). In practice however, what we can do is to sub- 
stitute the series expansion of a and b in terms of s and t given in (52) or to 
better accuracy in appendix B. This finally leads to the ultimate expansion 
of i<2 and i*3 in s and t which we seek. Let us present results after setting 

(64) T = 1 and z = t 2 . 

We can consider the perturbative expansion of F g in s with exact dependence 
on t or the perturbative expansion in t with exact dependence on s. In the 
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former case, we find 

(65) 6 

^M) = ^ + E^ 2 ), 

with the constant terms —1/240 and 1/1008 reproducing the free energies of 
the Gaussian model (48) in accordance with the fact that our model reduces 
to the Gaussian one for st = 0. 

We find that indeed C21 (z) = C%{z) and C$i(z) = C 3 (z) reproduce 
C g (z) given in (12), and C2,b(z) and Czj,{z), for b = 2, 3, reproduce gener- 
ating functions given in (18). This confirms that our model correctly repro- 
duces these known results for b = 1, 2, 3 backbones in genera two and three. 
Moreover, we can easily perform expansions to higher order in s revealing 
generating functions for more backbones. For example for four backbones, 
we find 

144z 7 

C 24 (z) = —.(38675 + 620648z + 2087808z 2 

(1 — iz) 6 

+ 1569328z 3 + 134208/), 
C 3A (z) = 48 ^i 6 (53416125 + 1194366915z + 6557325096z 2 
+ 10738411392z 3 + 4580024832z 4 + 236239616z 5 ), 
and for five backbones, we find 
144/ 

C 25 (z) = — - (2543625 + 62424520z + 375044396/ 

(1 — \z) iV l l 

+ 671666053/ + 314761848/ + 18335696/), 
720/° 

C 35 (z) = - (360380790 + 11275076865^ + 95744892585/ 

(1 - AzY 7 ' 2 

+ 282797424880z 3 + 291167707410/ + 85497242928/ + 3218434848/). 

We can also determine the perturbative expansion in t with exact depen- 
dence on s. Again scaling variables as in (64) and now expanding up to 
seven chords, i.e., to order z 7 , we find 

1 s 2 
F 2 = -— + 21s/ + (483s + 1485— J* 5 + (6468s + 26808s 2 + 15015s 3 )/ 

2 i 7Rec;ft/|„3 , OQOn^rio^ ^7 i /»V.v8\ 



+ (66066s + 526064s 2 + 768564s 3 + 232050s 4 )/ + 0{z 

1 s 2 

+ 1485s/ + (56628s + 225225— )z 7 + 0(/), 



1008 v 2 
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These expansions provide the numbers of configurations with a fixed number 
of chords and arbitrary numbers of backbones in accordance with Proposi- 
tion 1.1. For example, for genus two with seven chords, there are at most 
four backbones, and for genus three with seven chords, there are at most two 
backbones. We can easily continue the above expansions to higher powers 
of s and t. 



Appendix A. Genera zero and one free energies 

In this appendix, we determine free energies for genus zero and one, i.e., 
compute Fq and F\ in general in notation introduced in section 3, and 
ultimately substitute to derive results for our model in section 5. 

A.l. Genus zero. The genus zero free energy for the matrix model given 
in (20) can be written in a saddle point approximation [12, 29] as 

Fq = —T J dx p(x)V (x) + T 2 J dxdy p(x) p{y) log \x — y\ + 
+TA( J p{x)dx-t ) 

(66) = ~[j dxp(x)V(x)-Xt ). 

In this variational problem, it is crucial to impose the condition preserving 
the number of eigenvalues, which is fixed by (24). This can be done by 
introducing and subsequently determining the Lagrange multiplier A, where 
the rest of the notation in the above equation is the same as in section 3. 
The Lagrange multiplier can be determined from the asymptotic condition 
on y{x) which (for the one-cut solution) leads to the formula 

(67) A = lim (2 J y(x)dx - V(A) - i log A) , 

which is finite as singular terms from the upper integration limit must be 
canceled by V(A) and to log A in the limit of large A. Substituting the results 
found in section 5, we find 
(68) 

tD 2 + S{3t 2 D 2 + 2tS(5 - 3tS) - 4) + 2tD 2 {3tS - 1) log(-D/2) 

8t(tS - 1) ' 

and to defined in (24) takes form 

^ ^ (a - b) 2 (3(a + b)t - 2) _ D 2 (3tS - 1) 



16T((o + 6)t-2) 4T(tS-l) 



Moreover, we find the following integral 
(70) 

Jdxp(x)V(x) = - m s T ' ts _ 1) {fi + ^ts + tD-WS-tD-l) )' 



TOPOLOGICAL RECURSION FOR CHORD DIAGRAMS 



29 



where 

/i = 165 - 16(s + 3S 2 )t + 8S(-3D 2 + 2s + 6S 2 )t 2 + 

+ (D 4 - 16S 4 + AD 2 (2s + 7S 2 ))t 3 + 3D 2 S{D 2 - 4(2s + S 2 ))t 4 

h = 8(-25 + 2(s + AS 2 )t - 4S(-D 2 + s + 3S 2 )t 2 + 

+ 2{S 2 - D 2 ){s + AS 2 )t 3 + S(D 2 (3s + 4S 2 ) - 2D 4 - 2S 4 )t 4 ). 

Substituting now (68), (69) and (70) into (66), we find Fq expressed in terms 
of cut endpoints a and b or rather their combinations S and D defined in (53). 
We can now substitute expansions of S and D obtained from solving (51) 
or (54). In particular, substituting perturbative expansions of S and D in s 
with exact dependence on t as given in appendix B, setting T = 1 and z = t 2 
as in (64), and neglecting z-independent contributions which reproduce the 
genus zero result in the Gaussian case, we find perfect agreement with known 
results. 

Indeed, Cq,i(z) = Cq{z) reproduces the generating function for the Cata- 
lan numbers given in (10) while Co, 2(2) and Cq^(z) agree with (18). More- 
over, we can perform expansions in s to arbitrary order and find exact 
generating functions for arbitrary numbers of backbones. In particular for 
four and five backbones, we find 

_ . , 24z 3 (3 + 18z + 8z 2 ) 

Coa{z) ~ {T^Izy " 

(71) = 72z 3 + 2448z 4 + 44544z 5 + 585984z 6 + . . . 

24z 4 (55 + 741z + 1512z 2 + 336z 3 ) 
C °M ~ (l-4z)^ " 

(72) = 1320z 4 + 67944^ 5 + 1765440z 6 + . . . 

To test these results and to provide an idea of their combinatorial complexity, 
we show in appendix C by an explicit enumeration that in genus zero on 
four backbones there are indeed 72 diagrams with three chords and 2448 
diagrams with four chords. 

A. 2. Genus one. As shown in [14], the genus one free energy is given by 
(6), where M(x) is in general defined via (30) and in our model takes the 
form (56), where a and b are given by a solution to the system (51). We can 
expand this result to arbitrary accuracy in s or t to find generating func- 
tions for chord diagrams with arbitrary numbers of backbones or chords. 
The constant term is logarithmic in T and merely reproduces the Gaussian 
genus one free energy, Setting T = 1 in accordance with (64), expanding in 
s, keeping the exact dependence in t and using results presented in appendix 

B, we find that indeed C\^{z) = C\{z) agrees with C\(z) given in (12) and 
£1,2(2)) Ci^(z) agree with (18). This confirms that our approach correctly 
reproduces known results in genus one for three or fewer backbones. More- 
over, performing expansions to higher order in s, we can immediately find 
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generating functions for more backbones, for instance, 
24z 5 (715 + 7551z + 12456z 2 + 20962 3 ) 



(l-4z)!0 

48z 6 (13650 + 253625z + 975915z 2 + 840600z 3 + 86448z 4 ) 

(1 - 4z)25/2 

Appendix B. Cut endpoints 



It is convenient to present here a solution for cut endpoints a and b ob- 
tained from equations (51) in terms of variables introduced in (53) 

(73) S-2±», 

Equations (54) provide a solution which is perturbative in s and exact in t: 

_2TV2 _ 2st 2 T 1 / 2 3s 2 t 4 T 1 / 2 (3 + 4£ 2 T) 

(1 - 4t 2 r) 3 / 2 (1 - \t 2 Tf 

s 3 t 6 T i/2( 47 + 308t 2 T + 12 8t 4 T 2 ) 

(1 - 4i 2 T) 13 / 2 
s 4 t 8 T 1 / 2 (1051 + 15756i 2 T + 33168i 4 T 2 + 6720t 6 T 3 ) 
4(1 - 4t 2 T) 9 

3s 5 t 10 T 1 / 2 (2037 + 53448t 2 T + 262224t 4 T 2 + 272640t 6 T 3 + 32768i 8 T 4 ) 



4(1 - 4t 2 T) 23 / 



2 



+ 



st 2s 2 t 3 (l + 8t 2 T) 

~ (1 - 4t 2 T) 3 / 2 + (1 - 4i 2 T) 4 

s 3 t 5 (7 + \m 2 T + 280t 4 T 2 ) 
+ (l-4t 2 T) 13 /2 

2s 4 i 7 (15 + 684t 2 T + 3552t 4 T 2 + 2560i 6 T 3 ) 
(1 - U 2 Tf 

s 5 t 9 (143 + 10408t 2 T + 107094t 4 T 2 + 240144t 6 T 3 + 96096t 8 T 4 ) 6 
+ (1 - 4i 2 T)23/2 + 0{ - s > 

Appendix C. Explicit diagram counting 

To illustrate the combinatorial complexity of chord diagrams and fully 
appreciate the generating functions which we have computed, we directly 
count in this appendix the number of genus zero diagrams on four backbones 
with three or four chords. In this way, we also confirm correctness of the 
coefficients co,4(3) = 72 and co,4(4) = 2448 in the genus zero expansion on 
four backbones in (7). 

Our first task is to distribute endpoints of chords along backbones in all 
possible inequivalent ways. If there are n chords, then there are 2n such 
endpoints, perhaps better apprehended as half-chords. We shall denote a 



0{s t 
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partition of {1, 2, • • • 2n} into i%i sets of cardinality i > 1 as {ni, . . . , n^} 
where J2k=i n k = 2n. 

We must then count how many ways there are to produce a connected 
chord diagram of the correct genus by pairing up these putative endpoints. 
In fact in the cases under consideration, the genus is always zero by Propo- 
sition 1.1, so here we must count simply the number of connected such 
pairings. 

To begin, consider diagrams on four backbones with three chords. The 
partitions of six are: {3, 1, 1, 1} and {2, 2, 1, 1}. 

• For {3, 1, 1, 1}, all possible linearly ordered sets of backbones arise 
by embedding them in the real line with the induced linear order and 
orientations. To achieve such an embedding, we merely determine 
which among the four backbones has three half-chords. Order these 
half-chords 1, 2, 3, which connect uniquely to the other backbones 
occurring in the order 1, 2, 3 along the real axis. There are thus 4 
times 3! = 24 such chord diagrams. 

• For {2,2,1,1}, there are (2) = 6 places to put the two 2's. By 
connectivity, there must be a chord between the two 2's, and there 
are four ways to do this. The two remaining half-chords on the 2's 
again come in an order along the real axis as do the two remaining 
l's, giving another factor 2!, for a total of 6 • 4 • 2! = 48. 

The grand total is thus 24 + 48 = 72 diagrams in agreement with the first 
term in the expansion of the generating function Co,4 given in (7). 

We next consider the more involved case of four backbones with four 
chords. There are the five partitions {5,1,1,1}, {4,2,1,1}, {3,2,2,1}, 
{2, 2, 2, 2}, and {3, 3, 1, 1} of the eight half-chords into four backbones, and: 

• For {5, 1, 1, 1}, there are four locations for the 5. Two half-chords on 
5 must be connected, and there are (I) = 10 of these. The remaining 
half-chords on 5 connect arbitrarily to the l's giving another factor 
3!, for a total of 4 • 10 • 6 = 240 diagrams for {5, 1, 1, 1}. 

• For {4, 2, 1, 1} there are four locations for 4 and three for 2. There 
must be at least one 4-2 chord by connectivity and at most two by 
cardinality. 

If there is a unique chord 4-2, among the 4 • 2 such possibilities, 
then the other half-chord on the 2 must connect to a 1 by definition, 
and there are two of these; the remaining 1 then can connect to any of 
the three remaining half-chords; there are a total of 4-3-4-2-2-3 = 576 
in this sub-case. 

There are (2) = 6 times two possible pairs of 4-2 chords and two 
ways to produce a connected diagram connecting the remaining l's 
to the remaining half-chords on 4; there are a total of 4-3-6-2-2 = 288 
in this sub-case. 
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There are a total of 576 + 288 = 864 chord diagrams for {4, 2, 1, 1}. 

• For {3,2,2, 1}, there are four locations for 3 and three for 1. There 
must be at least one 3-2 chord by connectivity. 

If there is exactly one 3-2 chord among the 3 • 4 possibilities, then 
the other half-chord from the 2 must go to the other 2, with two 
possibilities. This gives 4 • 3 • 4 • 3 • 2 = 288 in this sub-case. 

If there are exactly two 3-2 chords, then their 2-endpoints must lie 
on different backbones by connectivity. There are thus (g) = 3 times 
eight such pairs of chords and a unique completion to a diagram with 
exactly two 3-2 chords. There are a total of 4 • 3 • 3 • 8 = 288 in this 
sub-case. 

There are 4-3-2 configurations of three 3-2 chords with a unique 
completion, so a total of 4 • 3 • 4 • 3 • 2 = 288 in the sub-case of exactly 
three 3-2 chords. 

There are a total of 288 + 288 + 288 = 864 diagrams for {3, 2, 2, 1}. 

• For {2, 2, 2, 2}, there can be no two chords with the same backbone 
endpoints by connectivity, so the first chord on the first backbone 
lands in one of six places and the second in one of four places on dif- 
ferent backbones. There are two connected completions of any such 
partial diagram, for a total of 6 • 4 • 2 = 48 diagrams for {2, 2, 2, 2}. 

• For {3,3,1,1}, there are ( 2 ) locations for the 3's. A fixed 1 must 
connect to a 3 in one of three possible ways times two for the two 
3's. Either this 3 connects to the other 1 in two possible ways and to 
the other 3 in three possible ways, or not, in which case, the other 
1 connects to the other 3 in three possible ways with two possible 
connected completions for a total of 6 • 6 • (2 • 3 + 3 • 2) = 432 diagrams 
for {3,3,1,1}. 

Overall, we find 240 + 864 + 864 + 48 + 432 = 2448 chord diagrams in 
agreement with the second term co,4(4) in the expansion of the generating 
function Co,4 given in (7). 
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